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Observing Quantum Tunneling in Perturbation Series 
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ABSTRACT 

We apply Borel resummation method to the conventional perturbation series of 
ground state energy in a metastable potential, V{x) = — gx^/A. We observe 
numerically that the discontinuity of Borel transform reproduces the imaginary 
part of energy eigenvalue, i.e., total decay width due to the quantum tunneling. 
The agreement with the exact numerical value is remarkable in the whole tunneling 
regime < g ^ 0.7. 



★ e-mail: hsuzuki@mito.ipc.ibaraki.ac.jp 
f e-mail: yasuta@mito.ipc.ibaraki.ac.jp 



The quantum tunneling is a purely non-perturbative phenomenon: This phrase 
has been widely accepted. Let us consider a simple quantum mechanical example 

Since the ground state is metastable in this potential, the eigenvalue is defined 
by the analytic continuation from g' < 0. Equivalently, the Schrodinger equation 

is defined with a rotated boundary condition, ijj^x) for x = e^'^^H, with t 
±cx) [1], if the "escape out" solution is taken. With this boundary condition, the 
semi-classical (WKB) calculation [1] gives the imaginary part of energy eigenvalue 
for g' < 1, 
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This is related to the total decay width F due to the quantum tunneling,''' F — 
—2lmE{g). The semi-classical result (2) vanishes to all order of the expansion 
on g, i.e., the tunneling effect is invisible in a simple expansion with respect to the 
coupling constant. 

On the other hand, the conventional Rayleigh-Schrodinger perturbation series 
of ground state energy is given by 

oo 

m-T.'ng'', (3) 

n=0 

where the first several coefficients read^ 

CO = 1/2, ci = -3/16, C2 = -21/128, C3 = -333/1024, • • • . (4) 

For g real, the simple truncated sum of the expansion (3) is of course real and does 
not produce the desired imaginary part. Therefore usually one needs some non- 



I We may characterize the tunnehng regime of the couphng constant space as the real part 
of energy is lower than the potential barrier height. The numerical solution shows < g ^ 
0.677 is the tunneling regime. 

§ We have prepared the exact perturbative coefficients c„ to n = 153, using the technique 
in [1]. 
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perturbative technique, such as the semi-classical approximation of Schrodinger 
equation or, of the Euclidean path integral, to estimate the tunneling amplitude. 
In the language of the Euclidean path integral, the exponent in (2) is given by 
the action of the bounce solution [2] and the 0{g) term is the two loop radiative 
correction around the bounce [3]. 

In this article, we re-examine the above common wisdom. We try to reproduce 
the imaginary part by solely relying on the conventional perturbation series (3). 
This kind of approach has been pursued by Kleinert and co-workers recently [4]. 
In their work, a variational perturbation method, which is equivalent to the order 
dependent mapping (ODM) method [5], was applied to resum the perturbation 
series. They showed numerically that an accurate value of the imaginary part was 
reproduced in the strong coupling regime g ^ 0.1. Later, the rigorous convergence 
proof of ODM was given [6] for > g^o — 0.1, where Qq is the convergence radius 
of the strong coupling expansion. 

ODM works quite well in the strong coupling regime. However it cannot take 
the place of the semi-classical approximation because it fails in the weak coupling 
regime g ^ 0.1. This is not satisfactory from the above motivation although the 
failure can be understood [6] as the effect of Bender-Wu singularity [1] on the 
higher Riemann sheet of the complex g plane. As the another undesirable feature 
of ODM, the answer is given as a root of a higher algebraic equation, hence does 
not allow a simple analytical characterization such as (2).^ 

We nevertheless do not believe the above disadvantage of ODM itself is the 
fundamental difficulty of an approach based on the perturbation series. In fact, 
from the extensive studies in the seventies [7], it has been known that the large 
order behavior of perturbation series is determined by the tunneling effect (at least 
in super- renormalizable cases) . Therefore it is rather natural to expect the leading 
semi-classical behavior (2) for <^ 1 can be extracted from the information of 

^ H.S. thanks T. Eguchi for his notice on this point: This was the original motivation of the 
present work. 
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perturbation series. The idea behind our approach is simple: If the large order 
behavior of perturbation series is determined by the lowest semi-classical contribu- 
tion (2), the lower order perturbation coefficients, which we can compute reliably, 
should contain the information on the higher order corrections. We shall see below 
this expectation is in fact correct. 

This possibility is important even practically because the systematic higher 
order correction to the lowest semi-classical calculation is not a simple matter, 
especially in many variable systems: One should include the correction due to 
multi-bounce type configuration, with an integration over the (quasi-)collective 
coordinates, and the perturbative expansion around the multi-bounce, and so on. 
Such a problem, even in quantum mechanics, is still under active current researches 
(see, for example [8]). On the other hand, the perturbative calculation is the best- 
established technique in quantum theory. 

Here we take a classical and conservative approach to this problem. Namely 

we apply to the perturbation series the Borel resummation method [7,9], which 
proceeds as follows: We first define Borel transform from the perturbative coeffi- 
cients Cn in (3), 



oo 

Cr, 



n=0 

Then the energy eigenvalue is defined by Borel integral 

oo 

E{g) = - f dze-''<^B{z). (6) 



However the imaginary part of energy eigenvalue implies the Borel transform de- 
velops a singularity on the positive real axis. The imaginary part (2) for g ^1 is 
reproduced by a fractional branch point 

B{z)^ f{zQ){zQ-zr + ---, (7) 

where f{zQ) — — 2-\/2/7r, zq — 4/3 and a — —1/2; this is the singularity nearest 
to the origin and there may exist other singularities. The integration along the 
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positive real axis (6) is therefore ill-defined as it stands (the so-called "non-Borel 
summable" case) and we define the integral by deforming the integration contour 
to the upper side of branch cut. From the final result, we shall see this choice of 
contour corresponds to our boundary condition in the original Schrodinger problem. 

Incidentally, to find the position of the nearest singularity of Borel transform zq, 
we have used the information of semi-classical calculation (2). In principle, those 
values of zq and a may be found solely from the perturbation series Cn- Assuming 
only the first term of (7), we have a relation 

{n + l)cn n+1 



Cn+1 n-a 

n + l)Cn nCn-1 



1-1 = 1 

ZQ I- {2a+l)/n+ a{a + l)/-n? 



whose n — > oo limit is nothing but the Appel's comparison theorem. We found 
numerically that zq — 1.33341 and a — —0.518130 for n — 152, which are more or 
less consistent with the exact values. After observing this, we will use the exact 
value zq =4/3 in what follows. 

What we have to do is a construction of Borel transform (5) from the perturba- 
tion coefficients Cn and Borel integration (6) with the deformed contour. Unfortu- 
nately this simple recipe does not work practically because the convergence radius 
of series (5) is finite (= zq) due to the branch point singularity (7). To perform 
the integration along the positive real axis, we have to continue analjd;ically the 
series (5) outside the convergence circle, which is impossible without having all the 
coefficients of the series. As is well known, however, the analytic continuation can 
be avoided by the conformal mapping technique [10]. 

We introduce a new variable A by 



A 



This transformation maps the whole cut z plane into an interior of a unit circle 
on A plane. In particular, an interval z e [0, zq) is mapped to A e [0, 1) and. 
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z E [zq + is, oo + ie] is to the upper arc of a circle with radius 1 — £ on A plane. 
The infinity ^ = oo is mapped to A = —1. In terms of A, the series (5) is expressed 
as 

Bw^f^A'. 4.i:(-i)--" /_'^|;' (4.0)°^. (10) 

k=0 n=0 ^ / \ / 

The point is that, assuming the absence of singularity of B(z) on the cut z plane* 
the convergence circle of the series (10) now is the unit circle, within which the 
whole cut z plane is mapped. Therefore we may use a truncated sum of the 
expansion (10) in the integration (6). This trick has been extensively used in 
practical applications of Borel resummation [7,9]. However, in our present case, 
one subtlety arises: We have to integrate (10) along its convergence circle, on which 
the (uniform) convergence of the series is not guaranteed in general. For us, this 
seems to be one reason why the Borel resummation has not been seriously applied 
to the tunneling phenomenon.^ Anyway we can now go on, hoping the convergence 
of the method. 

We then parameterize the unit circle as A = e*^. Prom (10), (9) and (6), we 
find that the A^th order approximation of the imaginary part is given by 

[Im E{g)U ="-^ [d9 exp (-'-^^—) J^rffcsinfc^, (11) 

'^-^ 

where is defined in (10). Note that this is solely expressed by the perturbative 
coefficients to A^th order (except the value of zq for which we may use the semi- 
classical method, or n ^ oo limit of (8)). Therefore this is the analytic expression 
of tunneling amplitude which we were looking for. In fact, a contribution from 

* The branch point of B{z) at z = zq is transformed to a single pole at A = 1. We investigated 
other possible singularities of B{z) on A plane by applying to the series (10) the Appel's 
comparison theorem (8) (c„ dn)- It does not converge to a definite value and indicates 
the other singularities on a positive real axis of z plane z> zq. We suspect this is an effect 
of the multi-bounce configuration. 

f The application of Borel resummation to (3) with <? < (anharmonic oscillator) is found 
in [11]. 

6 



an infinitesimal interval near 9 ~ in (11) is proportional to ~ exp(— 2:0/5'), the 
semi-classical behavior. Remaining integration may be regarded as its higher order 
corrections? 

In Fig. 1, we plotted the logarithm of the relative error of (11) to the exact 
imaginary part for several values of g. We observe a very rapid convergence and 
practically we may even use A?^ ~ 5 which only gives a few percent error for 
g ~ 0.3! In Fig. 2, eq. (11) with = 15 is plotted as a function of g. For a 
comparison, the exact numerical value and the semi-classical calculation including 
the two loop correction (2) are also plotted. In Table 1, several values in the 
weak coupling regime are listed. The agreement with the exact value in a whole 
tunneling regime < ^ 0.7 is remarkable and, as we expected, it is much better 
than the semi-classical result: It provides a systematical quantitative improvement 
of the bounce calculation. 

Wc have seen that for the mctastablc potential (1), the (generalization of) Borcl 
resummation method combined with the conformal mapping technique reproduces 
the very accurate tunneling amplitude. We expect a similar result holds for generic 
metastable potentials, for which the imaginary part is "physical." On the other 
hand, the tunneling phenomenon also exists in a bounded potential but with multi- 
minima (the quantum coherence). In this case also the perturbation series is non- 
Borel summable due to the tunneling [12] and the Borel resummation produces 
imaginary part. However since the potential is bounded, this cannot be regarded as 
a physical imaginary part and our approach, based on a conventional perturbation 
series around the "trivial vacuum," cannot directly be applied to such a bounded 
potential. The imaginary part will merely give the order of magnitude of barrier 

I It might look strange that we have a finite imaginary part from the integration on A plane, 
because an integration along the real axis — 1 < A < 1 (which does not give the imaginary 
part) may be added to close the integration contour, and we have assumed the absence of 
singularity inside the unit circle. In fact there is no contradiction because the change of 
variable (9) produces in the integrand a double pole at A = — 1 and a infinitesimal — 7r/2 
rotation aroimd A = — 1 can produce minus of the imaginary part. 

§ This unphysical imaginary part should be canceled by an imaginary part arising from the 
quasi-collective coordinate integration of multi-instantons [13]. 
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penetration, or typical time scale of an oscillation among different minima. 

Finally we comment on the generalization to higher dimensional models. The 
quantum mechanics (1) we have analyzed is of course so simple that everything can 
be computed numerically. However it is equivalent to a one-dimensional {D — 1) 
scalar field theory: 



We can expect a similar set of formulas as (10) and (11) (for the energy density) 
works for (12), if D < 4 or with an ultraviolet cutoff. For these cases, the model 
is super-renormalizable or ultraviolet finite and thus the renormalon [14], another 
known source of the Borel singularity, does not emerge. An analysis of the Eu- 
chdean bounce solution shows that the nearest singularity is, zq — 5.850448 for 
D = 2 and ZQ = 18.897251 for D = 3, with a = -D/2 [15]. What we have to 
compute is (appropriately renormalized) vacuum bubble diagrams as much higher 
orders as possible. Then an accurate bubble nuclcation rate in a certain D — 1 di- 
mensional system will be reproduced. A study along this line is in progress. 



The work of H.S. is supported in part by the Ministry of Education Grant-in- 
Aid for Scientific Research, Nos. 08240207, 08640347, 08640348 and 07304029. 
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Figure Captions 

Figure 1. The logarithm of the relative error of (11) to the exact value for, g — 0.08 
(the full circle), g — 0.3 (the circle) and g — 0.6 (the square), as a function 
of N. 

Figure 2. The imaginary part (normalized by the lowest semi-classical result, the first 
term of (2)) computed by (11) with iV = 15 (the full circle). For a compari- 
son, the exact numerical value (the solid line) and the semi-classical result (2) 
(the broken line) are plotted. To indicate the convergence behavior, the result 
with = 5 is also plotted for g < 0.16 by circles. 

Table Caption 

Table 1. Comparison of (11) [N = 15) and the exact numerical value of imaginary 
part in the weak coupling regime. All the numbers are normalized by the 
lowest semi-classical result, the first term of (2). 
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exact 


semi-classical (eq. (2)) 


eq. (11) with = 15 


0.08 


0.914996 


0.920833 


0.915012 


0.1 


0.891039 


0.901042 


0.891023 


0.12 


0.865072 


0.881250 


0.865082 



Table 1 
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